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LONG-TERM  GOALS 

Our  long-term  goal  is  to  understand  how  flows  in  near-coastal  zone  (20m  to  100m)  respond  to  a 
variety  of  forcing  mechanisms  including  tidal  pressure  gradients,  surface  waves,  surface  heating  and 
cooling,  surface  wave-bottom  current  interaction,  and  tidally  generated  bottom  boundary  currents. 
Because  the  nature  of  this  response  varies  throughout  the  water  column  and  depends  strongly  on  the 
non-linear  coupling  of  stratification,  turbulence  and  flow  structure  characterizing  the  structure  of  the 
water  column  in  this  environment  is  a  very  difficult  field  measurement  task. 


1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

30  SEP  2003 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2003  to  00-00-2003 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


4.  TITLE  AND  SUBTITLE 

A  Study  of  the  Structure  of  the  Near-Coastal  Zone  Water  Column  Using 
Numerical  Simulations 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Environmental  Fluid  Mechanics  Laboratory „Department  of  Civil  and  report  number 
Environmental  Engineering„Stanford  University„Stanford„CA,  94305 

9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Our  long-term  goal  is  to  understand  how  flows  in  near-coastal  zone  (20m  to  100m)  respond  to  a  variety  of 
forcing  mechanisms  including  tidal  pressure  gradients,  surface  waves,  surface  heating  and  cooling,  surface 
wave-bottom  current  interaction,  and  tidally  generated  bottom  boundary  currents.  Because  the  nature  of 
this  response  varies  throughout  the  water  column  and  depends  strongly  on  the  non-linear  coupling  of 
stratification,  turbulence  and  flow  structure  characterizing  the  structure  of  the  water  column  in  this 
environment  is  a  very  difficult  field  measurement  task. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

8 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


OBJECTIVES 


It  is  possible  to  gain  some  insight  into  the  physics,  and  into  our  ability  to  model  or  parameterize  the 
physics,  by  looking  at  a  more  idealized  version  of  this  problem  using  large  eddy  numerical  simulations 
(LES).  We  are  performing  simulations  of  the  mixing  processes  in  the  upper  layer  of  the  near-coastal, 
and  deeper,  using  a  periodic  channel  as  the  computational  domain.  Benefits  of  using  numerical 
simulations,  as  compared  to  laboratory  or  field  experiments,  is  the  relative  ease  with  which  information 
about  the  turbulence  can  be  extracted  from  the  flow  and  the  control  over  the  external  variables.  The 
study  has  two  fundamental  goals:  (1)  Developing  deeper  understanding  of  the  interaction  between 
various  physical  mechanisms  that  affect  the  dynamics  of  the  upper  mixed  layer  in  the  near-coastal 
ocean  and  the  deep  ocean;  and  (2)  Developing  improved  parameterizations  of  these  processes  for  use 
in  large  eddy  simulations  (LES)  and  Reynolds-Averaged  Navier-Stokes  (RANS)  models  to  be  used  in 
modeling  on  the  larger  scales. 

APPROACH 

The  LES  technique  is  employed  to  solve  the  Navier-Stokes  equations  for  stratified  turbulent  channel 
flow  numerically.  A  vortex  forcing  term  has  been  added  to  simulate  Langmuir  circulations.  The  code 
employs  the  dynamic  eddy  viscosity  model  of  Germano  et  al.  (1991)  with  Lilly's  (1992)  least  squares 
modification.  This  code,  which  was  originally  developed  and  implemented  on  the  400-node  Intel 
Paragon  XP/S  supercomputer  at  SDSC  by  Garg  et  al.  (1994),  has  been  modified  both  to  be  portable  to 
extant  parallel  computing  platforms  and  to  account  for  the  additional  physical  phenomena  we  wish  to 
examine.  The  code  was  run  with  a  128x256x128  grid  on  our  local  40-node  Compaq/Alpha  Beowulf 
cluster  (http://fluid.stanford.edu/baywulf/baywulf.html),  and  we  also  ran  the  code  with  similar 
resolution  on  the  260-node  IBM  Netfinity  Linux  Supercluster  at  the  Maui  High  Performance 
Computing  Center  (http://www.mhpcc.edu/doc/huinalu.html). 

WORK  COMPLETED 

Following  Zhou  (1999),  we  added  the  Craik-Leibovich  vortex  forcing  tenn  to  the  momentum 
equations  in  order  to  simulate  Langmuir  circulations.  We  also  added  the  ability  to  vary  the  thickness  of 
a  thermocline  density  profile  in  the  channel.  Due  to  increases  in  computing  power,  higher  Reynolds 
number  (ReT  =  480  and  720)  cases  have  also  been  completed  to  complement  the  Re  T  =  360  cases.  We 
also  modified  the  code  to  allow  us  to  run  either  Poiseulle  flow  cases,  Couette  flow  cases,  and  mixed 
Poiseulle-Couette  flows  with  and  without  Langmuir  forcing.  These  modifications  to  the  code  allow  us 
to  examine  the  interactions  among  Langmuir  circulations,  stratification,  and  smaller  scale  turbulence. 
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Finally,  we  used  the  DNS  simulations  reported  in  Holt  et  al  (1992)  and  Shih  et  al  (2000)  together  with 
laboratory  experimental  data  from  Barry  et  al,  (2001),  Rehmann  and  Koseff  (2003),  and  Jackson  and 
Rehmann  (2003)  to  examine  paramterization  of  the  eddy  diffusivity  in  stratified  turbulent  flows. 

RESULTS 

We  ran  numerical  simulations  of  stratified  channel  flows  for  a  number  of  different  flows,  varying  the 
strength  of  the  stratification,  the  type  of  forcing  (Poiseulle,  Couette,  or  Poiseulle-Couette),  the  relative 
position  of  the  stratification  in  the  water  column,  and  the  type  of  density  boundary  conditions  (either 
constant  temperature  or  zero  flux). 

Stratified  Channel  Flows:  We  were  able  to  compute  planar  averaged  statistics  for  these  flows  which 
are  typically  characterized  by  a  relatively  well-mixed,  lower  boundary  region  where  the  turbulence  is 
quite  vigorous;  a  thermocline  region  where  the  turbulence  is  strongly  affected  by  the  buoyancy  forces; 
and  an  upper  mixed  layer  which  is  separated  from  the  lower  boundary  by  the  thennocline.  The  more 
weakly  stratified  flows  (with  zero  flux  boundary  conditions)  are  quickly  mixed,  while  in  the  more 
strongly  stratified  cases  the  turbulence  is  significantly  suppressed  over  the  entire  water  column.  In  the 
free-stream  portion  of  the  flows,  in  general,  the  turbulence  behaves  similarly,  albeit  not  identically,  to 
turbulence  in  homogeneous  shear  flow.  The  nature  of  the  stratified  turbulence  is  very  similar  to  that 
seen  in  other  flows.  This  is  shown  in  Figure  1  where  the  density  and  turbulence  profiles  for  three 
different  flow  cases  (weak,  mild,  and  strong  stratification)  are  shown  as  a  function  of  time.  Rex=  360 
and  the  zero-flux  boundary  conditions  are  imposed. 
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Figure  1:  Time  evolution  of  density  and  turbulence  profiles  for  runs  with  Rer=  360,  and  a)  RiT=  1, 

b)  RiT  =  1,  and  c)  RiT  =60. 


Langmuir  Circulations:  We  have  completed  the  simulations  for  a  whole  range  of  flows  with 
Langmuir  circulations  using  the  stratified  channel  flows  (described  above)  as  the  basis.  To  generate  the 
Langmuir  cells  we  need  to  apply  an  additional  surface  stress  forcing.  These  results  are  currently  being 
analyzed  as  part  of  another  ONR  sponsored  project.  We  have  found  that  the  growth  of  Langmuir 
circulations  is  inhibited,  but  not  completely  suppressed,  by  stratification.  The  presence  of  the  vortex 
forcing  mechanism  above  a  thennocline  induces  surface  mixed  layer  deepening. 

Stratified  Shear  Flows:  We  have  used  the  results  from  our  very  extensive  database  of  homogeneous, 
sheared,  stratified  turbulence  to  look  at  parameterizations  of  the  eddy  diffusivity,  Kp.  W  found  that 
the  Osbom  (1980)  prediction  is  reasonable  when  the  turbulence  is  stationary  (20  <s/vN  <100)  but  for 
values  that  are  greater  than  100  the  diffusivity  appears  to  be  a  function  of  (s/vN  )  .  Furthermore, 
from  dimensional  analysis  the  diffusivity  requires  2  dimensionless  parameters  for  full  description:  the 
Froude  number  and  the  turbulent  Reynolds  number.  We  found  that  a  high  Reynolds  the  diffusivity 
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became  independent  of  Reynolds  number,  and  we  were  able  to  develop  a  relationship  for  Kp  in  terms 
of  Froude  number  alone. 

IMPACT/APPLICATIONS 

The  simulations  completed  demonstrate  the  intrinsic  value  of  DNS  and  LES  in  that  it  allows  us  to 
calculate  each  term  in  a  model  or  parameterization  of  the  extant  physics.  Evaluation  of  existing 
turbulence  closure  models  or  commonly  used  sub-grid-scale  parameterizations  is  therefore  a  lot  more 
complete  than  with  experiments  alone.  Our  simulations  of  the  channel  flows  are  the  first  important 
step  in  developing  a  code  for  studying  the  evolution  of  the  density  structure  of  the  water  column  in  the 
near-coastal  ocean.  Once  completed  this  code  will  be  a  valuable  research  tool  for  use  in  conjunction 
with  field  work  currently  underway  involved  in  measuring  flowfields  in  the  near-coastal  ocean. 

TRANSITIONS 

The  numerical  databases  developed  have  been  analyzed  by  the  Pi’s  in  other  research  projects  and  the 
data  has  been  used  by  researchers  at  other  institutions.  For  example  Diamessis  and  Nomura  at  UCSD 
have  used  the  data  from  the  simulations  of  Shih  et  al.  (2000)  to  further  examine  the  interaction 
between  vorticity  and  rate-of-strain  in  stratified  turbulence.  In  addition,  Ivey  at  the  Center  for  Water 
Research  at  the  University  of  Western  Australia,  Rehmann  at  the  University  of  Illinois,  and  Umlauf  at 
EPFL  in  Switzerland,  are  using  the  data  to  look  at  parameterizations  of  turbulent  length  scales  and 
buoyancy  flux  in  stratified  flows. 

RELATED  PROJECTS 

Shear  Production  and  Dissipation  in  a  stratified  tidal  flow  -  ONR  -  (Monismith  PI).  Our  field  work 
includes  work  on  stratified  tidal  flows  in  which  we  are  making  Reynolds  stress  measurements  using 
broad-band  ADCPs 

An  Experimental  Study  of  a  Breaking  Interfacial  Wave  -  NSF-  (Koseff  PI).  In  the  laboratory  we  are 
performing  experiments  in  an  attempt  to  measure  the  mixing  associated  with  a  breaking  internal  wave 
at  a  stratified  (two-layer)  interface  using  the  wave-generation  technique  of  Rapp  and  Melville.  In  this 
work  we  are  measuring  the  mixing  efficiency  associated  with  such  an  event. 
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